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^ ■ ABSTRACT 

| We present observations and CMB foreground analysis of the Parkes Galactic Meridian 

Survey (PGMS), an investigation of the Galactic latitude behaviour of the polarized 
synchrotron emission at 2.3 GHz with the Parkes Radio Telescope. The survey consists 
of a 5° wide strip along the Galactic meridian I = 254° extending from Galactic plane 
to South Galactic pole. We identify three zones distinguished by polarized emission 
' properties: the disc, the halo, and a transition region connecting them. The halo section 

qq , lies at latitudes |6| > 40° and has weak and smooth polarized emission mostly at large 

■ scale with steep angular power spectra of median slope /3 m od ~ —2.6. The disc region 

covers the latitudes \b\ < 20° and has a brighter, more complex emission dominated 
by the small scales with flatter spectra of median slope /3 mo d = —1-8. The transition 
region has steep spectra as in the halo, but the emission increases toward the Galactic 
| plane from halo to disc levels. The change of slope and emission structure at b ~ —20° 

is sudden, indicating a sharp disc-halo transition. The whole halo section is just one 
environment extended over 50° with very low emission which, once scaled to 70 GHz, 
\ is equivalent to the CMB i?-Mode emission for a tensor-to-scalar perturbation power 

j_j ■ ratio rhaio = (3.3±0.4) x 1CP 3 . Applying a conservative cleaning procedure, we estimate 

d ' an r detection limit of Sr ~ 2 x 10~ 3 at 70 GHz (3-sigma C.L.) and, assuming a dust 

polariztion fraction < 12%, Sr ~ 1 x 10~ 2 at 150 GHz. The 150 GHz limit matches 
the goals of planned sub-orbital experiments, which can therefore be conducted at 
this high frequency. The 70 GHz limit is close to the goal of proposed next generation 
space missions, which thus might not strictly require space-based platforms. 
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1 INTRODUCTION 

The study of the Galactic polarized synchrotron emission is 
essential for two cutting-edge fields of current astrophysics 
research: the detection of the B-Mode of the Cosmic Mi- 
crowave Background (CMB), for which the Galactic syn- 
chrotron is foreground emission, and the investigation of the 
magnetic field of the Galaxy. 

The CMB B-mode is a direct signature of the 



primordial gravitation al wave background (GW B ) left 
by inflation (e.g., Kamionkowski fc Kosowskvl 1 19981 : 



iBovle. Steinhardt, fc Turokl l2006t ). The amplitude of its 
angular power spectrum is proportional to the GWB power, 
which is conveniently expressed relative to the amplitude of 
density fluctuations — the so-called "tensor-to-scalar pertur- 
bation power ratio" r.Q Still unde tected, the current up per 
limit is set to r < 0.20 (95% C.L.. iKomatsu et alj|2009h by 
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1 Refer to lPeiris et~al] i2003l ) Eq. 10 for a full definition of r. 
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the results of the Wilkinson Microwave Anisotropy Probe 
(WMAP). A detection of the B-Mode would be evidence 
for primordial gravitational waves, and a measurement of r 
would help distinguish among several inflation models and 
investigate the physics of the early stages of the Universe. 

Reaching this spectacular science goal will be difficult 
because of the tiny size of the CMB B-Mode signal, fainter 
than the current upper limit of 0.1 fiK, and perhaps as faint 
as the 1 nK corre sponding to the smallest r acc essible by 
CMB (r ~ 10" 5 , lAmarie. Hirata fc Seliak1l2005l ). At such 
low levels the cosmic signal is easily obscured by the Galactic 
foreground of the synchrotron and dust emissions. 

Investigations of the synchrotron contribution have 
been conducted over recent years, but data are still insuf- 
ficient to give a comprehensive view (Fig. [TJ. IPaee et al.l 
(|2007h analysed the 23 GHz WMAP polarized maps and find 
that the typical emission at high Galactic latitude is strong: 
at 70 GHz it is equivalent t<0 r ~ 0.3, even higher than the 
cu rrent upper lim i t. An analysis of the same WMAP data 
bv lCarretti et all |2006b) identified regions covering about 
15 per cent of the sky with much lower emission levels, offer- 
ing a better chance for B-Mode detection. In these regions 
the polarized foreground is fainter, equivalent to a B-Mode 
signal corresponding to r in the range [1 x 10~ 3 , 1 x 10~ 2 ]. A 
better characerization of the polarized foreground is crucial 
especially for sub-orbital experiments (ground-based and 
balloon-borne), which will observe small sky areas. 

The design of future experiments is dependent on the 
frequency of minimum foreground emission. WMAP finds 
that th is is in the range 60-70 GHz for high Galactic lat- 
itudes dPage et aLl feoO?). but it remains unknown for the 
lowest emission part of the sky. It has been suggested that 
the dust emission might have deeper minima than the syn - 
chrotron in the areas of lowest emission (e.g.. lLangel f2008), 
shifting the best window for B-Mode detection to higher 
frequencies. 

Synchrotron emission from the Milky Way is not only 
a foreground for CMB polarization measurements, but can 
also be used to study the Galactic magnetic field. The total 
intensity of synchrotron emission can be used to estimate the 
total magnetic field strength, while the polarized intensity 
gives the strength of the regular component. This analysis in 
external galaxies has shown that the spiral arms are usually 
dominated by a small-scale, tangled, magnetic field with a 
weaker coherent large-scale field aligned with the arms. In 
the inter-arm regions the regular component dominates and 
in some spirals, magnetic arms with coherent scales up to 
the size of the disc have been detected in between the gas 
arms (e.g. see Beck 2008 for a review). 

The synchrotron emissivity of our own Galaxy is harder 
to understand because of our location inside it, but has the 
advantage that it can be studied in detail. Frequency de- 
pendent synchrotron depolarization can be used to deter- 
min e typical scale and str ength of small-scale magnetic fields 
fe.g lGaensler et ail [2001), and all-sky synchrotron emis- 
sivity maps can characterise the synchrotron scale height 
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Figure 1. Summary of the current knowledge of the .B-Mode 
power spectra of the Galactic synchrotron emission. Spectra are 
estimated at 70 GHz. A brightness temperature frequency spec- 
tral index of a = —3.1 has been assumed for all extrapolations. 
C*WMAP is the general contamination at high Galactic latitude, as 
estim ated by the WMAP team using a ~ 75% sky fraction at 22.8- 
GHz dPage et alj|2007h . Spectra measured in small areas selected 
for their low emission are also sho wn: the target fields o f the exper- 
iments BOOM ERanG (Cf^,, .ICarretti et alj|2005bj) and BaR- 
SPOrt fCgj . [Carretti et al.ll2006al ). and the up per limit found 
in th e two fields of the experiment DASI (C , p ASI ,[ 
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|200y). Finally, C.s/jy. --? shows the es timate of the emission in the 
best 15% of the sky jCarretti et alj|2006bl. the shaded area indi- 
cates the uncertainties), while Cdrao th e values found in som e 
intermediate regions using 1.4 GHz data jLa Porta et al.i r2006l. 
For comparison, CMB spectra for three values of r are also shown. 



l|Beuermann et aT1ll985h , or can be used for large-scale mod- 
eling of the Galactic magnetic field, especially in the halc0. 
The relative parity of the toroidal magnetic f i eld compo- 
nent i s still under discu s sion (see e.g.,lHanlll997l : iFrick et al.l 
l200ll : ISun et al.l[2008l ). Ijansson et aD (j2009r i use WMAP 
synchrotron maps at 23 GHz to show that the magnetic 
field behaviour in the Galactic disc and halo may differ con- 
siderably. 

Data from external galaxies does not help in constrain- 
ing the Milky Way magnetic halo, as there is a wide variety 
of magnetic field configurations: from galaxies without ev- 
ident halo field, to X-shaped fields centred at the galaxy 
centre, to large almost spherical magnetic halos (see Beck 
2008 for a review). 

Recent maps of pol arized Galactic synchrotron radia- 
tion at 1.4 and 22.8 GH z (|Wolleben et al.ll2006l ; lTestori et al] 
l2008l : IPage et al.ll2007h show polarized emission across the 
entire sky, and can be used to study the Galactic mag- 
netic field. However, the 1.4 GHz maps show that the disc 
emission is strongly depolarized up to lat itudes |6| ~ 30° 
jWolleben et al.l l2006l : iTestori et all 120081 ) . while Faraday 
depolarization effects are still present up to |&| w 40 — 50° 
ijCarretti et al.ll2005al ). Furthermore, those data consist of a 
single frequency band and do not enable rotation measure 



2 By equivalent to r we signify the strength of a foreground whose 
spectrum would match the spectrum of the CMB B-mode emis- 
sion at the £ =~ 90 peak arising from conditions characterised by 
the given value of r. 



3 Galactic halo is used in this context as the gaseous and mag- 
netic field distributions out of the Galactic disc, and is not nec- 
essarily connected to the stellar halo. 
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computations. The WMAP data at 22.8 GHz are virtually 
unaffected by Faraday rotation (FR) effects, but the sensi- 
tivity is not sufficient since, once binned in 2° pixels, about 
55% of the sky has a signal to noise ratio S/N < 3. This 
area corresponds to all the high Galactic latitudes with the 
exception of large local structures, which is most of the sky 
useful both for CMB studies and to investigate the Galactic 
magnetic field. 

Therefore, synchrotron maps at intermediate frequen- 
cies over all Galactic latitudes are needed to explore the 
behaviour of the contamination of the CMB with latitude 
as well as to study the Galactic magnetic field in the disc, 
the halo, and the disc-halo transition. 

In this work we present the Parkes Galactic Merid- 
ian Survey (PGMS), a survey conducted with the Parkes 
Radio Telescope to cover a strip along an entire south- 
ern Galactic meridian at 2.3 GHz. The area is free 
from large local structures, making it ideal for investi- 
gating both the CMB foregrounds and the Galactic mag- 
netic field. The PGMS overlaps the tar get area of severa l 
CMB e xperiments like BO OMERan G iMasi et al.| | 200rJl . 
QUaD jBrown et all 120091) BIC EP (|Chiang et all boOai , 
and EBEX (|Grainger et al.|[2008h . Our results may have di- 
rect implications for all these experiments. 

In this paper we present the survey, observations, and a 
characterisation of the polarized emission. We also present 
an analysis of the measured emission as a contaminating 
foreground to CMB B-mode studies. Analysis and impli- 
cations for the Galactic magnetic field will be subject of a 
forthcoming paper (Paper II, Haverkorn et al. 2010 in prepa- 
ration). A third paper will deal with the polarized extra- 
galactic sources (paper III, Bernardi et al. 2010 in prepara- 
tion). 

Survey and observations are presented in Section [2] the 
ground emission analysis in Section [31 and the maps in Sec- 
tion(4] The analysis of both the angular power spectrum and 
emission behaviour is presented in Section [S] while the dust 
contribution is investigated in Section [5] The detection lim- 
its of r are discussed in Section [7| and, finally, our summary 
and conclusions are reported in Section [SJ 



2 THE PARKES GALACTIC MERIDIAN 
SURVEY 

The available data and the properties of the synchrotron 
emission discussed in Section [T] lead to the following main 
requirements for a survey. Observations must 

(i) be conducted at a low enough radio frequency for the 
synchrotron emission to dominate the other diffuse emission 
components, but at a frequency higher than 1.4 GHz to avoid 
significant Faraday Rotation effects; 

(ii) cover all latitudes from the Galactic plane to the pole, 
to explore the behaviour with the Galactic latitude b; 

(iii) cover regions free from large local structures, such 
as the big radio loops, that would distort the estimates of 
typical conditions at high latitudes. 

The Parkes Galactic Meridian Survey (PGMS) is a 
project to survey the diffuse polarized emission along a 
Galactic meridian designed to satisfy these requirements. 



It surveys a 5° x 90° strip along the entire southern merid- 
ian I = 254° from the Galactic plane to the south Galactic 
pole (Fig. [2]). The observations have been made at 2.3 GHz 
with the Parkes Radio Telescope (NSW, Australia), a facil- 
ity operated by the ATNF - CSIRO Astronomy and Space 
Scienc^B a division of CSIRC0. It also includes an 10° x 10° 
extension centred at I — 251° and b = —35°. 

The selected meridian goes through one of the low emis- 
sion regions of the sky identified u sing the WMAP data 
(Fig. [2] see also lCarretti et al.ll2006bh and is free of large lo- 
cal emission structures. The meridian also goes through the 
area of deep polarization obse rvations of the BOOMERanG 
experiment (|Masi et al.l 1200(1); the 10° extension near b = 
—35° is positioned to best cover that field. 

At long wavelengths, measurements of Galactic polar- 
ized emission in regions of high rotation measure are cor- 
rupted by Faraday depolarization. At 1.4 GHz Faraday de- 
polar ization is significant u p to Galactic latitudes |6| < 40- 
50° (|Carretti et all l2005al ) where RM > 20 rad/m 2 . At 
2.3 GHz this RM limit increases to 60 rad/m 2 , allowing 
a clear view of polarized emission over all high Galactic lat- 
itudes and well into the upper part of the disc. 

The observations were made in four sessions from Jan- 
uary 2006 to September 2007 with the Parkes S-band Galileo 
receiver, named after NASA's Jupiter exploration probe for 
which the Parkes te lescope and this rece iver were used for 
down-link support (| Thomas et al.l Il997h . The receiver re- 
sponds to left- and right-handed circular polarization, whose 
cross-correla tion gives Stokes parameters Q and U (e.g. 
lKraus|[l986l ). This scheme provides more protection against 
total-power (gain) fluctuations than the alternative: a re- 
ceiver responding directly to the linearly polarized signals. 

The original feed used for the Galileo mission has been 
replaced by a wide-band corrugated horn, highly tapered to 
reduce sidelobes and the response to ground emission. The 
feed illuminates the dish with a 20 dB edge taper and the 
first side-lobe is 30 dB below the main beam. 

The ATNF's Digital Filter Bank 1 (DFB1) was used to 
produce all four Stokes parameters, /, Q, U , and V . DFB1 
was equipped with an 8-bit ADC and configured to give 
256 MHz spectrum with 128 2-MHz channels. Spectra are 
generated using polyphase filters that provided high spectral 
channel isolation. The isolation between adjacent channels 
is 72 dB, an enormous improvement over the 13 dB isola- 
tion of Fourier-based correlators. This, in combination with 
the high sample precision, gives excellent protection against 
RFI leaking from its intrinsic frequency to other parts of the 
measured spectrum. This is valuable in the 13-cm band as 
strong RFI can be present Q Recording spectra with spectral 
resolution greater than required for the polarimetry analysis 
has allowed efficient removal of RFI-effected channels, max- 
imising the effective useful bandwidth. Data were reduced 
to 30 8-MHz channels. The RFI removal typically yielded 
an effective total bandwidth of 160 MHz. 

The source B1934-638 was use d for flux c alibra tion as- 
suming the polynomial model by iRevnoldi (| 19941 ') for an 



4 http://www.atnf.csiro.au 

5 http://www.csiro.au 

6 The first observing session, in September 2005, used a Fourier- 
based correlator, and spectra were strongly contaminated by RFI. 
Subsequent use of DFB1 greatly improved the measurements. 
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Figure 2. The PGMS strip (dark gray) plotted on the WMAP polarized intensity map (L = yj ' Q 2 + U 2 ) at 22.8 GHz l lPage et a,l.ll2007t) 
binned in 2° pixels (HEALPIX pixelation with Nside=32). Pixels with S/N < 3 have been blanked (white). The PGMS goes through 
a region clear of large emission structures avoiding contamination by major large local anomalies, such as the big radio loops and the 
fan region (centred at I ~ 135°, b ~ 0°). The big radio loops are structures large several tens of degrees, like loop-I which extends from 
the north Galactic pole down to the Galactic plane, with a possible continuation into the southern hemisphere. The map is in Galactic 
coordinates with longitude I = 0° at centre and increasing leftward. 



accuracy of 5%. The polarization response was calibrated 
using the sources 3C 138 and PKS 0637-752, whose polar- 
ization states were determined using the Australia Telescope 
Compact Array (ATCA) with an absolute error of 1°. The 
statistical error of our polarization angle calibration is 0.5°. 
The astronomical IAU convention for polarization angles is 
used: angles are measured from the local northern meridian, 
increasing towards the east. It is worth noting that this dif- 
fers from the convention used in the WMAP data, for which 
the polarization angle increases westwards. The unpolarized 
source B1934-638 was also used to measure the polarization 
leakage, for which we measure a value of 0.4%. The off-axis 
instrumental polarization due to the optics response is about 
1%. 

The use of a system with both Q and U as corre- 
lated outputs mitigates gain fluctuation effects. To check the 
level of a 1/f noise component in the data we observed the 
South Celestial Pole (SCP), thereby avoiding azimuth (AZ) 
and elevation (EL) dependent variations of atmospheric and 
ground emissions. Remaining variations in the signal arise 
from intrinsic atmospheric changes and receiver fluctuations. 
Power spectra of the Q and U time-series are almost flat with 
no evidence of a 1/f component down to 3 mHz (see Fig. [3j| . 
This confirms that the system is stable and characterised 



by white noise up to 7-min time scales, sufficient for the 
duration of our scans. 

The Galactic meridian was observed in 16 5° x 5° fields 
and one 10° x 10° field. The fields are named PGMS-XX, 
where XX is the Galactic latitude of the field centre. Each 
field except PGMS-02 includes a 1° extension along b at the 
north edge for an actual size of 5° x 6° and an overlap of 
5° x 1° with the next northern field. 

The fields were observed with sets of orthogonal scans 
to give I- and 6-maps (scans along Galactic longitude and 
latitude, respectively). Each field was observed with 101 
latitude scans (fa-maps) and 121 longitude scans (Z-maps) 
spaced by 3 arcmin to ensure full Nyquist sampling of the 
beam (FWHM = 8.9 arcmin). The same sample spacing was 
used along each scan by scanning the telescope at 3°/min 
with a 1-second integration time. 

One full set of I— and fo-maps were observed for the 
6 disc fields at latitude \b\ < 30° and the 10° x 10° field, 
giving final a sensitivity of ~ 0.5 mK per beam-sized pixel. 
Over the ten high latitude fields (|6| > 40°), where a weaker 
signal was expected, two full passes were made to give a 
sensitivity of ~ 0.3 mK per beam-sized pixel. 

Prior to map-making, a linear baseline fit was removed 
from each scan and the ground emission contribution was 
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3 GROUND EMISSION 
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Figure 3. Power spectra of Q (top) and U (bottom) time series 
for a south Celestical pole observation. Both spectra are mostly 
flat with no evidence of a 1/f component up to 3 mHz. 



Table 1. Main features of the PGMS observations. 



Central frequency 
Effective bandwidth 
Useful bandwidth 1 
FWHM 

Channel bandwidth 
Central Meridian 
Latitude coverage 
Area size 
Pixel size 
Observation runs 



Q, U beam-size pixel rms sensitivity (halo fields) 
Q, U beam-size pixel rms sensitivity (disc fields) 



2300 MHz 
240 MHz 
160 MHz 
8.9 arcmin 

8 MHz 
l = 254° 
b = [-90°, 0°] 
5° x 90° 
3' x 3' 
Jan 2006 
Sep 2006 
Jan 2007 
Sep 2007 
0.3 mK 
0.5 mK 



After RFI channel flagging. 



estimated and cleaned up by the procedure described in Sec- 
tion [31 

The map-making procedure is based on the algorithm 
by Emerson & Grave (1988), which combines I— and fo-maps 
in Fourier space and recovers the power along the direction 
orthogonal to the scan, otherwise lost through the baseline 
removal. The algorithm is highly efficient and effectively re- 
moves residual stripes. 

Table [T] summarises the main features of the PGMS 
observations. 



Ground emission can seriously affect continuum observa- 
tions, especially in our halo fields where the emission has 
a brightness of only a few mK. Our tests have shown that 
the highest ground emission occur in the Zenith cap at eleva- 
tions above 60° where large fluctuations are observed in the 
data. Even though not yet fully understood, the most likely 
reason is the loss of ground shielding by the upper rim of the 
dish at large elevations. The receiver is located at the prime 
focus and is shielded by the dish up to this elevation, receiv- 
ing only atmospheric contributions from beyond the upper 
rim. Above this limit the ground becomes visible, contribut- 
ing a ground component that rapidly varies as the telescope 
scans. To avoid this contamination all PGMS observations 
were limited to the elevation range EL = [30°, 60°], between 
the lower limit of the telescope's motion and this region of 
ground sensitivity. 

Even though these precautions have significantly re- 
duced the effect, some contamination is still present in the 
halo fields, requiring us to develop a procedure to estimate 
and clean the ground contribution. 



3.1 Estimate and cleaning procedure 

The procedure is based on making a map of the ground 
emission in the AZ-EL reference frame. Any AZ-EL bin 
gathers the contributions of data taken at different Galactic 
coordinates and the weak sky emission is efficiently averaged 
out. Since the PGMS meridian goes through low emission 
regions, its high latitude data are ideal for such an aim. 
The smooth behaviour of the ground emission enables the 
use of large bins, which further helps average out the sky 
component. We therefore use a bin size of AEL = ldeg in 
EL and average over 8 degrees in Azimuth. The binning 
is performed in the instrument reference frame before the 
correction for parallactic angle. 

For a given 8 degree AZ bin we smooth the map along 
the EL direction. This reduces the residual local devia- 
tions that are mainly due to strong point sources. We use 
a quadratic running fit: for each EL bin, we fit the 7 bins 
centred at it with a 2-degree polynomial, and the bin value 
is then replaced with the fit result at the bin position. 

In the Azimuthal direction the data are sufficiently 
smooth that no fit is necessary in that dimension. We there- 
fore shift our 8 degree azimuth averages by 1 degree incre- 
ments, performing the elevation fit for each 1 degree bin. 
This results in a map of the ground emission in the AZ- 
EL frame with a bin size of 1 degree in both Azimuth and 
Elevation. 

The ground emission contamination mainly comes from 
the far lobes, which are frequency dependent. We generate 
ground emission maps for each frequency channel. A set of 
these maps is generated at each observing session. They are 
checked for constancy of ground emission before a grand 
average set is formed from all observations. 

Fig. U shows example AZ cuts in the 2300 MHz spec- 
tral channel. Over a range of about 30° in both AZ and 
EL the ground emission varies by less than 50 mK, smaller 
by an order of magnitude than ground emi ssion variations 
repor ted for other polarization surveys (e.g. Wollcbcn ct al. 
120061 ). We attribute this low response to ground emission to 
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Figure 4. Examples of azimuth cuts of the ground emission maps 
of both Q (top) and U (bottom) from the 2300 MHz channel of 
the September 2007 observations. 



the high edge tapering of the S-band feed and its consequent 
small sidelobes. 

Finally, the sky emission measurements, observed on 
a 3-arcminute grid, were cleaned of ground emission using 
the model just described. For each sky measurement, the 
ground emission at its actual (AZ, EL) was obtained by 
linear interpolation through a stan dard Cloud-In- Cell tech- 
nique (|Hocknev fc Eastwoodlll98l1 ). 



3.2 Cleaning procedure tests 

The low emission fields in the halo have been observed twice 
giving two independent maps taken in different runs and at 
different AZ and EL. That way, they are contaminated by 
different ground emission enabling us to test the cleaning 
procedure and estimate the residual contamination. The er- 
ror map can be estimated as half the difference of the two 
maps, in which the sky is cancelled, leaving the noise and 
any residual (unmodelled) ground emission. 

An example is given in Fig. [5] in which one of the fields 
with lowest emission is shown: PGMS-87. The images show 
Stokes Q (left) and U (right) of both the observed (top) and 
difference maps (bottom). 

The difference maps are clearly dominated by white 
noise, indicating that most of the ground emission has been 



removed at the level needed to measure the sky signal, and 
that the residual ground emission does not contribute sig- 
nificantly to the error budget. That residual can be seen as 
faint shadows of large scale structure in the difference maps. 

These visual impressions of the ground emission removal 
we know quantify by measuring the angular power spectr43 
of both the sky and the difference map. Fig. [S] reports the 
mean of the E— and _B-mode power spectra - (E+B)/2 
- which is the most complete description of the polarized 
emission. 

The spectrum of the sky signal is dominated by the 
diffuse emission at low multipole £, where it follows a power 
law Ce oc £ ,B with a steep slope (/? < —2.0). A flattening 
occurs at the high-£ end due to both noise and a point source 
contribution. A white noise spectrum would be flat (Ce = 
constant). 

The difference spectrum also follows a power law, but 
is much flatter than the sky signal. The best fit slope is 
/Snoiso = —0.70 ±0.06, which, although not pure white noise, 
is close to the ideal /3 no i SB — 0. Furthermore, the difference 
between sky signal and noise increases at large angular scale, 
giving a rapidly increasing S/N. 

The ground emission has a smooth behaviour that 
makes the largest scales the most susceptible to contami- 
nation. On a scale of 2° the rms fluctuation of the difference 
map is N20 — 60 /iK, a factor 2.5 larger than that expected 
from white noise (24 /iK), but much smaller than the sky 
signal (few mK). 

This field (PGMS-87) is the worst with regard to ground 
emission residuals; over all fields, angular spectral slopes 
/3noise fall in the range [-0.7, 0.0] and the rms noise on 
the 2-degree scale A^2° lie in [24, 60] /iK. The mean val- 
ues over all fields are /3 no i sc = —0.4 and N20 — 43 /iK. Once 
the pure white noise component is subtracted, the effective 
contribution by only ground emission can be estimated in 
^Vgrnd,2° = 36 /iK. With such results the impact of the 
ground emission may be considered marginal on the final 
mapping. 



4 PGMS MAPS 

Maps of the Stokes parameters Q, U and the polarized inten- 
sity L of the PGMS meridian are shown in Fig. andO 
while Fig. \W\ displays the whole 10° x 10° field (PGMS- 
34). All images are smoothed with a Gaussian filter of 
FWHM = 6' to give a better idea of the sensitivity on beam- 
size scale, for an effective resolution of FWHM = 10.7'. All 
data at latitude |6| > 30° are plotted with the same in- 
tensity range to show clearly the power and morphological 
differences. The disc fields (\b\ < 30°) require a more ex- 
tended scale. The two strongest sources present in our data 
(Pic A and NGC 612 in field PGMS-34 and PGMS- 77, re- 
spectively) have been blanked before the map generation 
of their fields. Without blanking, the high brightness range 
causes the map-making procedure to generate artefacts. 

The disc has the strongest emission, extending to lat- 
itude |6| ~ 20° with little variation of emission power. At 



7 See Section \E\ for the description of the power spectrum com- 
putation. 
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Figure 5. Q (left) and U images (right) of the observed (top) and difference maps (bottom) of the field PGMS-87 (6=[-84°, -90°]). 
No smoothing is applied for a resolution of FWHM=8.9'. Position coordinates are Galactic latitude and longitude; the brightness unit is 
mK, the intensity colour scales are linear. 
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Figure 6. Power spectrum of the mean (E+B)/2 of E— and 
B— Modes measured for the PGMS-87 polarized emission: sky 
emission (diamonds), difference map (stars) and best fit curves 
(dashed). 



higher latitudes, the emission starts to decrease up to the 
halo where it settles on levels one order of magnitude lower. 
The clear visibility of the bright polarized disc emission 



and its contrast with the fainter halo is a new result, not 
apparent from previous observations carried out at lower 
frequencies where the disc is strongly depolarized up to 
\b\ ~ 30°. This allows us to locate the disc-halo boundary 
in polarization, which has not been visible so far because of 
either strong depolarization (at 1.4 GHz) or insufficient sen- 
sitivity (at 23 GHz). A more quantitative analysis is given 
in Section [5] but the visual inspection of the maps clearly 
shows that it starts at |6| ~ 20°. 

The emission of the halo has a smooth behaviour with 
the power mostly residing on large angular scales. The disc 
emission is also smooth, at least at latitudes higher than 
|&| = 6°-7°. Closer to the Galactic plane the pattern has a 
more patchy appearance, suggestive of Faraday depolariza- 
tion effects being significant at 2.3 GHz. 

This supports the view that Faraday depolarization ef- 
fects are marginally significant in the disc, and are relevant 
only in a narrow belt a few degrees wide across the Galactic 
plane. 

Several polarized point sources are visible, especially in 
the halo where the diffuse emission fluctuations are smaller. 
To enable a cleaner analysis of the diffuse component the 
sources are identified, fitted and subtracted from the maps. 
Each source is located by a 2D-Gaussian fit of the stronger 
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Figure 7. Q (left), U (mid), and polarized intensity L = y Q 2 + U 2 images (right) of the six PGMS fields in the latitude range 
b = [-30°, 0°] (PGMS-27 through PGMS-02). Position coordinates are Galactic latitude; the brightness unit is mK. 
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Figure 8. As for Fig. [7] but in the range b = [-60°, -30°] (PGMS-57 through PGMS-34: of the latter only the 5° across the meridian 
I = 254° are imaged). Position coordinates are Galactic latitude; the brightness unit is mK. 
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Figure 9. As for Fig. but in the range b = [-90°, -60°] (PGMS-87 through PGMS-62). Position coordinates are Galactic latitude 
and longitude; the brightness unit is mK. 
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Figure 10. Q (top), U (mid), and L images (bottom) of the whole 10° X 10° fields (PGMS-34). Position coordinates are Galactic 
latitude and longitude; the brightness unit is mK. 
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Table 2. Polarization flux limits Sp lm used to selected the polar- 
ized sources. The chosen limits are latitude dependent. 



6-rangc 


S« m [mjy] 


[-10°, 0°] 


40 


[-20°, -10°] 


30 


[-30°, -20°] 


20 


[-40°, -30°] 


15 


[-90°, -40°] 


10 



component, either Q or U. Its position is then fixed in the fit 
of the second and weaker component to improve the fit ro- 
bustness. A polarization flux limited selection is applied with 
threshold set to ensure S/N ratios of at least 5. The ampli- 
tude of fluctuations in the maps is dominated by sky emis- 
sion (rather than by the instrument sensitivity), which varies 
along the PGMS meridian. The threshold we use is therefore 
a function of Galactic latitude, running from 10 mjy at high 
latitudes up to 40 mjy near the Galactic plane (Table |2J. 

In this work, the point source identification is carried 
out only for cleaning purposes. The catalogue and a detailed 
analysis of their properties are subject of a forthcoming pa- 
per (Bernardi et al. 2010, in preparation). 



5 ANGULAR POWER SPECTRA 

The angular power spectra (APS) of E— and B-Mode of 
the polarized emission have been computed for each field. 
They account for the 2-spin tensor nature of the polarization 
and give a full description of the polarized signal and its 
behaviour across the range of angular scales. In addition, 
the E— and B-Modes are the quantities predicted by the 
cosmological models enabling a direct comparison with the 
CMB. 

To cope with both the non-square geometry and the 
blanked pixels at the locations of the two brightest sources, 
we use a method based on the two-point correlation func- 
tions of the Stok es parameters Q and U described by 
Sb arra et all <|2003h . The correlation functions are estimated 
on the Q and U maps of the regions as 



C x {8) =XiXj 



X = 



■ U 



(1) 



where X% is the emission in pixel i of map X, and i and 
j identify pixel pairs at distance 8. Data are binned with 



pixel-size resolution. Power spectra C. 
integration 



are obtained by 



C 



C 



wi 



\C q {9)F 1 m{8) + C u {6)F 2 , e2 {6)] sm(9)d9 (2) 



W[ I [C u (0)F 1 , e2 {S)+C Q (e)F 2 ,i 2 (S)}^(G)de (3) 



where - Fj.im and F 2 j m a re functions of Legendre polynomi- 
als fsee lZaldarriagalll998 j for their definition), and Wf is the 
pixel window function accounting for pixel smearing effects. 

Since the emission power is best described by the quan- 
tity £(£ + l)C^/(27r), hereafter we will denote an angular 
spectrum following a power law behaviour Ct oc I® as 



• flat, if P = 
angular scales; 



-2.0: power equally distributed across the 



• steep, if P < —2.0: large scales dominate the power bud- 
get; 

• inverted, if /3 > —2.0: small scales dominates. 

We tested the procedure using simulated maps gener- 
ated from a known input power spectrum b y the proce- 
dure synfast of the software package HEALPix ( Gors ki et al.l 
2005). The input spectra are power laws with different 
slopes; for each slope we generated 100 simulated maps and 
compute their APS. The mean of the 100 APS reproduced 
the input spectrum and its slope correctly, with the excep- 
tion of an excess at the largest scales, mainly at the first 
two multipole bands. E— and B-Mode are related to the 
polarization angle pattern and this excess is likely due to 
the discontinuity of the pattern abruptly interrupted at the 
area borders. To account for this we corrected our spectra 
for the fractional excess estimated from the simulation as the 
ratio between the mean of the computed and input spectra. 

For a cleaner measure of the diffuse component, the 
point sources are subtracted from the polarization maps. 

The E— and B-mode spectra Cf and Cf have 
been computed for the 17 fields along with their mean 
c (e+b)/ 2 _ + cf)/2. Artificial fluctuations are gen- 
erated on E and B spectra because of the limited sky cover- 
age of the individual areas, but their mean suffers less from 
that effect and is a more accurate estimator if the power is 
distributed equally between the two modes, as is the case 
for Galactic emission. In addition, C\ gives a full de- 

scription of the polarized emission which the two individual 
spectra cannot give separately. Therefore we mostly use the 
mean spectrum (E+B)/2 to investigate emission behaviour 
and properties in the following analysis. 



Fig. M shows c[ E+B) ' 2 for all the fields. As an ex- 
ample of all three spectra, Fig. [T^] shows those of the two 
fields PGMS-52, which is from the low emission halo, and 
PGMS-34, our biggest field and the area observed by the 
BOOMERanG experiment. All spectra are shown without 
correction for the window functions. In most fields the spec- 
tra follow a power law behaviour that flattens at the high 
multipole end because of the noise contribution. Exceptions 
are the four fields closest to the Galactic plane (PGMS-17 
to PGMS-02) where a power law modulated by the beam 
window function dominates everywhere. . 

We fit the angular power spectra to a power law modu- 
lated by the beam window function Wf for the synchrotron 
component and a constant term N for the noise: 



( -'2O0 



(—) 

1200/ 



W e B + N 



W 



(4) 



where C 200 is the spectrum at £ = 200 and X — 
E, B, (E+B) /2 denoting E-Mode, B-Mode, and their mean 
(E+B)/2. Possible residual contributions by point sources 
are accounted for by the constant term. 

Plots of the best fits are shown in Fig. Illl and ll2l while 
the parameters of the synchrotron component are reported 
in Table El 

To analyse the behaviour of the synchrotron compo- 
nent, the power law component of the best fit spectra are 
plotted together in Fig. [15] where they are also extrapo- 
lated to 70 GHz for comparison with the CMB signal. To 
determine the spectral index for the frequency extrapola- 
tion we computed a map of spectral index of the polarized 
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Figure 11. Angular power spectra Cj, of the 17 PGMS fields. Both the measured spectra (diamonds) and the best fit cv. 

(solid) arc plotted. 
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Figure 12. Angular power spectrum of E-Mode (top), B-Mode (mid), and their mean (E+B)/2 of the two fields PGMS-52 (left) and 
PGMS-34 (right). Both the measured spectra (diamonds) and the best fit curve (solid) are plotted. 



Table 3. Best fit amplitude C^, (referenced to i = 200) and angular spectral slope /3 X of the PGMS fields (X = E, B, (E + B)/2 
denoting E-, B-Mode, and (E+B)/2, respectively). 



Field 




Gfoo [^ k2 1 




Qfoo [MK 2 ] 


/3 s 




< +s)/2 [^] 


p(E+B)/i 


PGMS 


02 


4500 ± 280 


— 1.69 ± 0.04 


4860 ± 880 


-1.90 ±0 


13 


4610 ± 540 


-1.76 ±0.09 


PGMS- 


07 


3850 ± 380 


— 1.77 ± 0.07 


4120 ± 340 


-1.75 ± 


06 


4030 ± 340 


-1.76 ± 0.06 


PGMS- 


12 


2800 ± 200 


-1.86 ±0.05 


3160 ±410 


-1.82 ± 


09 


3070 ± 240 


-1.85 ± 0.06 


PGMS- 


17 


2090 ± 260 


-2.00 ±0.10 


1270 ± 140 


-1.67 ±0 


08 


1840 ± 190 


-1.91 ± 0.08 


PGMS- 


22 


550 ± 100 


-2.54 ±0.17 


246 ± 39 


-2.60 ±0 


16 


376 ± 31 


-2.39 ± 0.08 


PGMS- 


27 


287 ± 58 


-2.41 ±0.24 


131 ± 19 


-2.17 ±0 


19 


263 ± 20 


-2.51 ± 0.09 


PGMS- 


34 


183.5 ±8.5 


-2.54 ±0.07 


171.8 ± 9.0 


-2.46 ±0 


07 


180.0 ±7.1 


-2.51 ± 0.06 


PGMS- 


42 


79 ± 11 


-2.76 ±0.19 


61.9 ±6.3 


-2.24 ±0 


14 


74.1 ±4.0 


-2.57 ±0.08 


PGMS- 


47 


59.8 ±7.7 


-2.91 ±0.20 


49.7 ±4.3 


-2.56 ± 


13 


56.3 ±3.1 


-2.78 ± 0.08 


PGMS- 


52 


47.9 ±4.4 


-2.74 ±0.13 


42.7 ±5.9 


-2.36 ± 


17 


51.3 ±3.2 


-2.66 ± 0.08 


PGMS- 


57 


46.7 ±3.4 


-2.62 ±0.10 


16.6 ±2.2 


-2.54 ±0 


21 


32.4 ± 1.7 


-2.61 ± 0.07 


PGMS- 


62 


28.9 ±3.1 


-3.27 ±0.22 


19.1 ±3.4 


-3.35 ± 


40 


24.1 ± 1.9 


-3.27 ±0.17 


PGMS- 


67 


15.6 ±2.0 


-2.17 ±0.19 


5.4 ± 1.1 


-2.23 ± 


34 


12.6 ± 1.6 


-2.35 ± 0.20 


PGMS- 


72 


40.1 ±4.1 


-3.08 ±0.16 


14.9 ±2.8 


-2.81 ± 


33 


29.7 ±2.7 


-3.03 ± 0.15 


PGMS- 


77 


41.6 ±6.4 


-3.05 ±0.27 


33.8 ±3.6 


-2.47 ±0 


16 


38.4 ±4.5 


-2.75 ± 0.19 


PGMS- 


82 


16.2 ±3.3 


-2.36 ±0.35 


18.0 ±3.0 


-2.51 ± 


27 


21.7 ±2.8 


-2.57 ±0.21 


PGMS 


87 


26.9 ± 1.6 


-2.21 ±0.07 


11.3 ±0.7 


-1.74 ±0 


07 


19.9 ± 1.4 


-2.08 ±0.09 
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Figure 13. Map of frequency spectral index of the polarized syn- 
chrotron emis sion computed using t he all-sky polarizati on maps 
at 1.4 GHz iWolleben et al.l |2006|; iTestori et al.l |2008|) and at 
22.8 GHz CWMAP-5vr. Iffinshaw et al.ll2009l 'l 
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Figure 14. Distribution of the spectral indexes reported in 
Fig. [13] 



synchrotron emission (Fig. I13[) using the all - sky polariza- 
tion s urveys at 1.4 GHz ( Wolleben et al.ll2006l; iTestori et alj 
I2008T ) and 22.8 GHz (WMAP. iHinshaw et alj|2009h . The in- 
dex distribution at high Galactic latitudes (|6| > 30°) peaks 
at Qfpoi = —3.21 with dispersion Aa po i = 0.15 (Fig. I14[) . 
T his is consistent w ith the analysis of the WMAP-5yr data 
bv lGold et ail |2009() . who find a polarized synchrotron spec- 
tral index of Qwmap = —3.1 i n the WMAP frequency range, 
and that of iBernardi et al.l {2004) of total intensity data 
(ai — —3.1 in the range 1.4-22.8 GHz). We therefore as- 
sume a spectral index of Qf syn ch = —3.1 for extrapolations 
up to the CMB frequency window. This assumption is some- 
what conservative and gives some margin to our conclusions. 

Distinguishing fields in terms of their spectral slopes, 
Fig. [15] and Fig. [16] show the presence of two well defined 
regions: 

• The high and mid latitude fields (b = [-90°, -20°]), with 
steep spectra (f3 < —2.0); slopes are distributed over the 
wide range /3 ~ [-3.0, -2.0] (except for a few outliers). There 
is no clear trend with latitude and the slopes are rather 
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Figure 15. Best fits to the spectra of i?-Mode (top), B— Mode 
(mid), and their mean (E+B)/2 (bottom) of all PGMS fields. 
The plot reports only the power law component which describes 
the synchrotron emission. All PGMS fields are plotted together 
for a direct comparison of the behaviour with the Galactic lati- 
tude. The colour code goes from blue throughout red by a rain- 
bow palette for the areas from the south Galactic pole (PGMS- 
87) throughout the Galactic plane (PGMS-02), respectively. The 
spectra are scaled to 70 GHz for a comparison with the CMB 
signal (frequency spectral index a = —3.1). CMB B-Mode spec- 
tra for different values of tensor-to-scalar power ratio r are shown 
for comparison. The quantity plotted here is t{t+ l)/(27r) * , 
which provides a direct estimate of the power distribution through 
the angular scales: a flat spectrum means the power is evenly 
distributed; a decreasing spectrum (steep) is dominated by the 
largest scales; a rising spectrum (inverted) is dominated by the 
smallest scales. 
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Table 4. Mean (j3 x ), median (3* ,), and dispersion (erg) of the 
angular slopes at mid-high Galactic latitudes for all the three 
spectra X = E, B, (E + B)/2. 
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Table 5. As for Table[5]but for low Galactic latitudes. 
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Figure 16. Angular spectral slopes of the PGMS fields plot- 
ted against the field's centre latitude for i?-Mode (top), B— Mode 
(mid), and (E+B)/2. 



uniformly distributed. The mediarjf|is P^ B)/2 = -2.6 (Ta- 
ble[SJ). The dispersion a /3 ( E +B-)/2 ~ 0.24 is significantly larger 
than the individual measurement errors, meaning that this 
wide spread is an intrinsic property of the synchrotron emis- 
sion at these latitudes. 



We prefer to use the median to estimate the typical angular 
slope in this region because of the possible significant deviations 
by the outliers. 



• The low latitude fields (b = [-20°, 0°]), which show 
inverted spectra (/3 > —2.0); the slopes lie in a much nar- 
rower range ( / g( B + s '/ 2 = [—1.90, —1.75]) with dispersion 
o" /3 (e+b)/2 = 0.08 (Table [SJ. All spectra have mean and me- 
dian slope fi x = /3 x cd — —1.8, which can be considered the 
typical value of this region. 

This change from steep to inverted spectra is quite sud- 
den and clearly separates two different environments: the 
mid-high latitudes, characterised by a smooth emission with 
most of the power on large angular scales, and the disc fields, 
whose power is more evenly distributed with a slight pre- 
dominance of the small scales. 

Does this change indicate an intrinsic feature of the po- 
larized emission of the disc, or is it the effect of Faraday 
modulation, which transfers power from large to small an- 
gular scales in the low latitude fields? The answer is un- 
clear with the information available, but some points can be 
noted. In the disc the ISM is more turbulent than in the halo 
and the intrinsic emission might have more power on small 
angular scales. In addition, the low-latitude lines of sight go 
through much more ISM including more distant structures; 
these are expected to give more power to the small angular 
scales. Also the smooth emission of the two highest latitude 
disc fields (PGMS-12 and PGMS-17) make the presence of 
significant Faraday depolarization unlikely. Finally, we have 
computed the power spectrum of the individual frequency 
channels to search for a possible variation of the angular 
slope with frequency. Since the lowest frequencies would be 
more affected, steeper spectra at highest frequencies would 
support the presence of FR effects. We find that all the four 
disc fields have non-significant slope variation compatible 
with zero within 1.0-1.5 sigma, with the only exception of 
PGMS-02 which approaches 2-sigma. All these points sup- 
port the view that the structure of the low-latitude polarized 
emission derives from the intrinsic nature of the synchrotron 
emitting regions close to the plane, and is not imposed by 
Faraday depolarization along the propagation path. 

Considering the amplitude distribution of the PGMS 
fields, we further divide the mid-high latitude region iden- 
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Figure 17. As bottom panel of Fig. 1151 but with spectra 
grouped according to latitude region. The disc fields (solid, red: 
b = [— 20°, 0°]) are clearly distinguished by their spectral slopes 
and higher amplitudes, and the fields in the transition region (dot- 
ted, black: b = [—40°, —20°]) have amplitudes quite distinct from 
both the disc and halo fields (solid, blue: b = [-90°, -40°]). 



tified above into an halo (b — [—90°, —40°]) and transition 
region (6 — [—40°, —20°]). Thus we identify three distinct 
latitude sections: two main regions (disc and high latitudes) 
well distinguished by both emission power and structure of 
the emission, and an extended transition about 20° wide 
connecting them. Fig. [17] reports the spectra for all field, 
showing how they belong to the three regions, which: 

• Halo (b = [—90°, —40°]): the emission is weak here and, 
scaled to 70 GHz, is between the peaks of CMB models 
with r = 10 -3 and r = 10" 2 . The weakest fields (PGMS-87 
and PGMS-67) match models with r = 10~ 3 . The fluctua- 
tions from field-to-field dominate with no clear trend with 
latitude. A weak trend might be present with the emission 
power increasing toward lower latitudes, but the effect is a 
minor in comparison to the dominant field-to-field fluctua- 
tions. 

• Galactic disc (6 = [— 20°, 0°]): the emission is stronger, 
about two orders of magnitude brighter than that of the 
halo. Within the area there is no large variation of the emis- 
sion power, but slight increase toward the Galactic plane is 
evident. 

• Transition strip (b — [—40°, —20°]): here a transition 
between the faint high latitudes and the bright disc occurs. 
This is clearer at large scales, where the northernmost field 
(PGMS-22) is almost as bright as the weakest disc field 
(PGMS-f7) and the southermost field approaches the up- 
per end of the halo brightness range. 

An important consequence is the identification of a clear 
transition between disc and halo. The sudden change in the 
angular spectral slope at \b\ ~ 20° and the approximately 
constant emission power from the Galactic plane up to that 
transition clearly separate the 20° equatorial zone from the 
higher latitudes. Characterised by a more complex structure 
of the ISM this area can be associated with the Galactic disc. 

A second environment characterised by steep spectra 
and low emission is clearly present for \b\ > 40°. Both an- 
gular slope and amplitude exhibit wide fluctuations without 
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Figure 18. Mean of the power spectra of the halo fields (b = 
[-40°, -90°]) for £-Mode (top), B-Mode (mid), and their mean 
(E+B)/2 (bottom). Both the mean of the measured spectra (di- 
amonds) and its best fit curve (solid) are plotted. 

any clear trend with latitude. We consider this high Galactic 
latitude section as a single environment with regard to its 
polarized synchrotron emission properties. Characterised by 
a smoother emission and simpler ISM structure, this area 
we associate with the Galactic halo. 

The emission of the halo section is very weak. In spite 
of large fluctuations, once scaled to 70 GHz the synchrotron 
component is equivalent to r values between 1 x 10~ 3 and 
7 x 10~ 3 , which matches the weakest areas observed so far in 
polarization. It is worth noticing that PGMS fields PGMS- 
87 and PGMS-67 have the weakest polarized synchrotron 
emission observed so far. 

The high Galactic latitudes above 40° are thus just one 
environment, at least in a low emission region not contami- 
nated by local anomalies like the area intersected by PGMS. 
This is very important for CMB investigations, since it tells 
that, in principle, it is possible to find large areas with op- 
timal conditions (extended over 50°, in the PGMS case). 

It is thus important to measure the mean emission prop- 
erties of the entire halo section (|b| > 40°). The mean spec- 
trum of the 10 halo fields and its best fit are plotted in 
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Table 6. Best fit amplitude and angular spectral slope /3 X of the mean spectrum of the PGMS halo section (X = E, B, (E + B)/2 
denotes E-, S-Mode, and (E+B)/2, respectively). 



Field Cj m [ M K 2 ] d E Cf 00 [^] pB C^ +B) ' 2 [ M K 2 ] 0(b+b)/2 

PGMS HALO 40.7 ±3.2 -2.65 ±0.10 28.5 ± 3.0 -2.64 ±0.12 35.5 ± 2.9 -2.60 ± 0.09 



B — Mode PGMS halo Dust: Temperature 




( E -\- B ) /2 

Figure 19. Best fit of the mean halo spectrum scaled 
to 70 GHz (solid) alongside CMB B— Mode spectra for different 
r values (dashed). The mean synchrotron contamination at all 
high Galactic latitudes in also shown for comparison (C^ MAp ) 
together with the previous estimates in other low emission re- 
gions: Cg /N<2 (shaded area), C D rao> CfoOM' and ^BaR ( see 
Fig. [15] for details). 



Figure 20. Temperature power spectrum of the dust emis- 
sion model at 94 GHz in the halo section of the PGMS (b = 
[— 90° , — 40 deg]). The spectrum is the mean of those computed 
in the individual PGMS fields. Both the mean of the computed 
spectra (diamonds) and its best fit curve (solid) are plotted. The 
latter is a power law with slope /3j = —2.50 ± 0.14 (see text and 
Table 0. 



Fig. 1181 Table [6] gives the best fit parameters; the extrapo- 
lation to 70 GHz is shown in Fig. [T5] 

The angular slope is f) ~ —2.6 for all the three spec- 
tra E, B, and (E+B)/2 and is thus be considered the typ- 
ical slope of the halo section. Note that this matches the 
slope measured a t high latitudes by WMAP at 22.8 GHz 
jPage et alj|2007h . indicating that the power distribution 
through the angular scales is the same at 2.3 and 22.8 GHz. 
This further argues against the significance of Faraday de- 
polarization, which would have transferred power from large 
to small scales and made the angular spectra frequency de- 
pendent. 

Once scaled to 70 GHz, the amplitude is equivalent to 

r h aio = (3.3 ±0.4) x 10" 3 , (5) 

roughly in the middle of the range covered by the individual 
fields. As mentioned earlier, this is a very low level and corre- 
sponds to the faintest areas observed previously, which thus 
seem to be more the normal condition of the low emission 
regions rather than lucky exceptions. 

Finally, the field PGMS-34 deserves mention as the area 
observed by the experiment BOOMERanG. It lies in the 
transition region, and although at the high-latitude end of 
this zone, has emission about five times greater than fields 
in the halo. It is suitable for measuring the stronger i?-mode 
(the aim of the 2003 BOOMERanG flight), but our results 
identify more suitable fields for detecting the B-mode. 



6 DUST EMISSION IN THE PGMS HALO 
SECTION 

The dust emission is the other most significant foreground 
for CMB observations. It has a positive frequency spectral 
index and dominates the foreground budget at high fre- 
quency. 

An estimate of the local dust contribution in the same 
portion of halo covered by the PGMS is thus important to 
understand the overall limits of a CMB B-Mode detection 
in that area. 

However, no polarized dust emission has been detected 
over the PGMS region, and even the total intensity dust 
map of WM AP is noise dominated in that area. We there- 
fore use the iFinkbeiner. Davis fc Schlegell (|l999l ) model of 
the total intensity dust emission applying an assumed po- 
larization fraction. First, we generate maps at 94 GHz using 
their model-8 for each of the 10 PGMS fields at \b\ > 40°. 
The temperature power spectra of each are then computed 
and averaged together to estimate the mean conditions of 
the whole 5° x 50° section. Fig. [20] shows both the mean 
spectrum and its best fit, whose parameters are reported in 
Table [7] Within the errors, the angular slope matches well 
that of the synchrotron. 

The total polarized spectrum is estimated from this 
temperature spectrum assuming a polarization fraction 
,/poi = 0.10, as inferred for high Galactic latitudes 
from the Archeops experiment data l|Benoit et al.l 120041 : 
iPonthieu et al.l [20051 ). The spectrum is further divided by 
two to account for an even distribution of power between 
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Table 7. Best fit amplitude C^qq and angular spectral slope /3j 
of the mean Temperature dust spectrum at 94 GHz in the halo 
section of the PGMS (|6| > 40°). 



Field 


c? 00 [mk 2 ] 




PGMS halo Dust 


(12.0 ± 1.4) x icr 6 


-2.50 ±0.14 
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Figure 21. Estimates of the dust polarized emission in the halo 
section of the PGMS at 70 GHz (solid). The best fit to the tem- 
perature spectrum is used assuming a dust polarization fraction 
/pol = 0.10 and then scaled to 70 GHz. The synchrotron emission 
spectrum estimated in Section [5] is shown for comparison (dot- 
ted), as well as the CMB B— Mode spectra for three different r 
values. 



E— and B-Mode, a reasonable assumption for the Galactic 
emission. 

For frequency extrapolations, we use a single Planck 
function modulated by a power law Td oc v ad /(ekT — 1) 
with index ad — 2.67 and tempe r ature T = 16.2 K, which 
reproduces the iFinkbeiner et all (| 19991 ) model well in the 
range [70, 150] GHz. It is worth noticing that at 94 GHz 
that function is well approximated by a power law with in- 
dex aJ; L = 1.5, consistent with the 5-yr WMAP result of 
p = 1.8 ± 0.3 ± 0.2 (statistical and systematic error, 



iGold et al.ll2009l ). 

Our estimate of the B-Mode polarized dust spectrum at 
70 GHz is given in Fig. 1211 At this frequency the two com- 
ponents of Galactic polarized emission are approximately 
equal, and so the total polarized foregound is at a mini- 
mum. This frequency is mildly dependent on assumed po- 
larized fraction of the dust component: alternate assump- 
tions of five percent or twenty percent shift the frequency of 
minimum to 80 GHz and 60 GHz respectively. 

This result is similar to that of the general high Galactic 
latitudes and suggests that the synchrotron-to-dust power 
ratio is only marginally dependent on the strength of the 
Galactic emission. 



7 LIMITS ON r 

To estimate the detection limits of r in the presence of the 
foreground contamination in the PGMS halo section, we 



consider an experiment with resolution # cm b = 1° at the 
CMB frequency channel to have adequate sensitivity at the 
I ~ 90 peak. 

We also account for the cleaning provided by foreground 
se paration technique s. We consider the two cases discussed 
bv lTucci et ail l|2005l ): 

(i) cleaning by subtracting the foreground map scaled to 
higher frequencies using just one frequency spectral index 
for all pixels. It is a coarse method and represents a worst 
case. The residual contamination depends on the spread of 
spectral indices in the area. We use Aft, = 0.15 for the 
synchrotron componen t (Section [5] see also iBernardi et al.l 
12004 IGold et aljhoogh . and assume the same dispersion for 
the dust emission (Actd = 0.15). We refer to this method as 
clean #1. 

(ii) cleaning by foreground subtraction, but assuming 
knowledge of the frequency spectral index for any individual 
pixel. In this case the amplitude of the residual contamina- 
tion depends on the measurement error of the frequency 
slopes. Here we assume a combination of the PGMS data 
with a synchrotron channel at 22 GHz onboard the CMB 
experiment with resolution S — 1° and sufficient sensitivity 
to give S/N = 5. For the dust emission, we assume that the 
CMB experiment includes a dust channel at 350 GHz with 
a sensitivity to allow a S/N = 5 and resolution scaled from 
the CMB channel (#d = flcmb^cmb/^d, where ^ cm b and v& 
are the frequencies of the CMB and dust channel) . We refer 
to this method as clean #2. 



These two methods are at the two ends of cleaning capabil- 
ities (#1 is coarser and less efficient, #2 is finer and more 
efficient), so give a good idea of the range of possible per- 
formance. 

The left panel of Fig. 122 1 shows the residual contamina- 
tions left by these two meth ods at 70 GHz in t he PGMS halo 
section using the method of lTucci et all (|2005h . The effect of 
method #1 is to reduce the amplitude of the contamination 
but preserve the shape of the power spectrum (for instance, 
the synchrotron residual preserves the angular spectral in- 
dex of -2.6). Method #2 gives comparable results at the 
I ~ 90 CMB peak, but has a flat white-noise-like spectrum 
which performs much better on large angular scales. While 
the simpler clean #1 looks appropriate as the target is the 
peak at I ~ 90, the clean #2 is better suited for the reion- 
ization bump at larger scales. 

With such low levels of residual foreground contribu- 
tion the effects of gravitational lensing become a dominant 
part of the contamination budget. Its subtraction is there- 
fore also required to realise the benefit of the low Galactic 
emission of the PGMS strip. Gravitational lensing effects can 
be cleaned using high resolution data (10 arcmin or better, 



ISeliak fc Hiratal 120041) and here we assume that it can be 
reduced by 10 fold jSeliak fc Hiratal liool ). 

The Fisher inform ation matrix is used to es t imate 
the detection limits (|Tegmark. Taylor fc Heavens! Il997l : 
Tcgmark 1997). If r is the only parameter to be measured, 
(other cosmological parameters being provided by Temper- 
ature and _E-Mode spectrum from other experiments such 
as WMAP or PLANCK), the Fisher matrix reduces to the 
scalar 
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B-Mode: cleaned (70 GHz) B-Mode: cleaned (150 GHz) 




Figure 22. Residual contaminations after foreground cleaning in the PGMS halo section at 70 (left) and 150 GHz (right) for synchrotron 
(red) and dust emission (blue). Results for both the two cleaning methods described in text are shown. The gravitational lensing 
contribution is also shown assuming that cleaning has reduced its contribution by a factor 10 in the power spectrum (green). 



Table 8. Detection limits of r at 3-sigma C.L. (Sr) at 70 and 150 GHz in the PGMS halo area and in a 2500-deg 2 region assumed to 
have equivalent foreground contamination (see text). Results for both cleaning types described in the text are reported. 



Area 


clean type 


Sr (70 GHz) 


Sr (150 GHz) 


PGMS 


clean #1 


2.2 X 10~ 3 


8.0 x 10~ 3 


PGMS 


clean #2 


1.5 x 10~ 3 


5.5 x 10~ 3 


2500-deg 2 


clean #1 


1.2 x 10~ 3 


4.2 x 10" 3 


2500-deg 2 


clean #2 


0.5 x 10~ 3 


1.8 x 10~ 3 



B\2 



dr 



and the rms error on r is 



(6) 



(T) 



The uncertainty ACf of the _B-mode spectrum is a function 
of the CMB spectrum Cf' clnb and the cleaning residuals of 
synchrotron, dust, and gravitational lensing: 



AC/ 



(2£ + l)A£f sky 



C 



C, 



B. synch — res 



i s~iB . dust — res . B , lens — res \ / \ 

+ O e +U e I , (8) 



where A£ is the width of the multipole bins and / s k y is 
the sky coverage fraction. As a set of cosmological param- 
eters we use the best fits of the 5-yr WMAP data release 
l|Komatsu et alj|2009h . 

At 70 GHz the detection limits of r in the PGMS halo 
region is 5r ~ 2 x HP 3 (3 -sigma C.L.) for both cleaning 
methods (Table [SJ. This is a very low level which makes the 
PGMS strip an excellent target for CMB experiments and 
enables accessing levels of r much lower t han previously es- 
timated for areas of comparable size (e.g. iTucci et al.ll2005l ; 
IVerde et alj I2006L who used higher foreground levels esti- 
mated from total intensity data). An important point is that 
there is only a marginal benefit from using the more sophis- 
ticated cleaning method. For a 250-deg 2 area most of the 
sensitivity resides in the I ~ 90 peak, where the dominant 



residual is the gravitational lensing and a better cleaning of 
the other contaminants is not critical. 

This result is therefore quite robust since it is based on 
actual measurements of the foreground contamination in a 
specific area and is marginally dependent on the cleaning 
method. Moreover, the leading residual term is the gravi- 
tational lensing, giving a good margin against errors in our 
dust polarization fraction assumption. 

It is also important to estimate Sr for 150 GHz, a fre- 
quency that, although far from the foreground minimum, 
is preferred by experiments based on bolometric detectors. 
However, as shown in Fig. 1221 right panel, the major residual 
at this frequency is dust emission, making the result depen- 
dent on the assumed dust polarization fraction / po i- A value 
of / P oi = 0.10 gives the limit Sr = 6-8 x 10" 3 (Table©, 
while the goal of the most advanced sub-orbital experiments 
planned for the next years (r~lx 10 -2 ) is still achievable 
under the reasonable assumption that / po i < 0.12. 

The PGMS halo region is a narrow strip 50° long and it 
is unlikely that its optimal conditions are confined to its 5° 
width. Therefore we consider that a larger area, say of 50° x 
50° extent, could be identified with properties comparable 
to those of the PGMS halo region. Such an area, about 6 
per cent of the sky, matches in size the southern portion 
of the low emission r egion identified in the WMAP data 
l|Carretti et al.ll2006bl ). 

The detection limit achievable over such an area drops 
to Sr = 5 x 10" 4 (3 -sigma C.L.) if method #2 is applied, or 
Sr = 1.2 x 10" 3 under the coarser method #1 (Table [8}. 
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Note that in this case there is a significant difference 
between the two cleaning methods, justifying the use of the 
more sophisticated method #2. This is because the larger 
area makes measurement of the lowest multipole compo- 
nents relevant, where method #1 is not effective. 



8 SUMMARY AND CONCLUSIONS 

The PGMS has mapped the radio polarized emission at all 
Galactic latitudes in a 5 degree strip at a frequency suffi- 
ciently high not to be affected by Faraday depolarization 
and with sufficient sensitivity to detect the signal in low 
emission regions. It is the largest area observed so far at 
high Galactic latitude uncontaminated by large local struc- 
tures. 

This has enabled us to investigate the behaviour of the 
polarized emission with latitude by computing the polarized 
angular power spectrum in 17 fields from the Galactic plane 
to the South Galactic pole. We can distinguish three lati- 
tude sections: two main regions well distinguished by both 
brightness and structure of the emission (disc and halo), and 
an extended transition connecting them. In detail they are: 

1. The halo at high Galactic latitudes (b = [-90°, -40°]) 
characterised by low emission fields with steep spectra (an- 
gular slope p = [—3.0, —2.0]), that is smooth emission 
dominated by large scale structures. The slope is almost 
uniformly distributed within a wide range, with median 
/3med = —2.6. 

2. A transition region at mid-latitudes (b — [—40°, —20°]) 
whose angular spectra are steep like those of the halo, but 
shows an increase of the emission power with decreasing 
latitude. 

3. The disc at low latitudes (b = [—20°, 0°]) characterized 
by inverted spectra with slopes in a narrow range of median 
/Smcd = —1.8. The amplitudes are two orders of magnitude 
brighter than in the halo and the power gradually increases 
towards the Galactic plane. 

The change in the angular slope around b = —20° is 
abrupt and identifies a sharp disc-halo transition from the 
smooth emission of the mid-high latitudes to the more com- 
plex behaviour of the disc; this is likely related to the more 
turbulent and complex structure of the ISM in the disc. 

The halo section has no clear trend with latitude of 
either emission power or angular slope, and, at least along 
the meridian sampled by PGMS, can be considered a single 
environment. This is very important for CMB investigations, 
as it indicates that it is possible to find large areas with 
optimal conditions for seeking the B-Mode. 

The synchrotron emission of the whole halo section is 
very weak. Once scaled to 70 GHz it is equivalent to r = 
3.3 x 10~ 3 , so that an experiment aiming for a detection limit 
of Sr = 0.01-0.02 would need no synchrotorn foreground 
cleaning. 

The dust component is also faint and equal to the syn- 
chrotron emission at 60-80 GHz for polarization fractions 
between 5 and 20 %. The frequency of minimum foreground 
in this low emission region is thus similar to that found with 
WMAP for the general high Galactic latitudes (75% of the 
sky). If confirmed in other regions, this would imply both 
that the dust-synchrotron power ratio is rather independent 



of the brightness of the Galactic emission and the frequency 
of minimum foreground nearly independent of the sky posi- 
tion. 

We estimate the r detection limit of this area account- 
ing for the use of foreground cleaning procedures. We apply 
both a coarse and a more refined method. The Galactic emis- 
sion is so low that the dominant residual contamination is 
from gravitational lensing, even assuming a 10-fold reduc- 
tion of the lensing foreground from cleaning. For both the 
two cleaning methods the detection limit is Sr ~ 2 x 10 -3 
(3-sigma C.L.) if the CMB B-Mode search is conducted at 
70 GHz. This result provides a sound basis for investigat- 
ing the B-Mode. The detection limit we have found here is 
even better than the goals of the most advanced sub-orbital 
experiments (r ~ 0.01, e.g.: SPIDER EBEX, and QU IET, 
ICrill et al.|[2008l : lGrainger et al.ll2008l : ISamtlebeifeOOSft and 
proves that there exists at least one area of the sky where it 
is realistic to carry out investigations of the B-mode down 
to very low limits of r. 

At 150 GHz the detection limit rises, but is still bet- 
ter than Sr = 0.01 assuming a reasonable dust polarization 
fraction (/ po i < 12%). 

These results are valid in the area actually observed by 
our survey. However, the PGMS halo section is extended 
over 50° along one dimension and it is unlikely that its ex- 
cellent conditions are confined to its 5° width. We have ex- 
plored the likely results from a larger 50° x 50° area, having 
properties similar to those of the PGMS halo section. In such 
a region the r detection limit would drop to Sr — 5 x 10~ 4 
(3-sigma) at 70 GHz. 

It is worth noticing that the gravitational lensing needs 
to be cleaned to take advantage of the low Galactic emission 
of the PGMS halo section. This can be effect ively carried out 
only with high resolution data (10' or finer. ISeliak fc Hiratal 
l2004h and the design of CMB experiments should comply 
with that rather than be limited to 1-degree to fit the peak 
at £ ~ 90. 

The results obtained here might suggest a review of 
plans to detecting the CMB B-Mode and associated inves- 
tigations of the inflationary scenarios. While the detection 
limit is limited to Sr = 1 x 10 -2 by the detector array size 
and sensitivity, observations at 150 GHz might be sufficient 
if conducted in an area like the PGMS with clear advantages 
of using the currently best detectors (bolometers) and of an 
experiment more compact than at 70 GHz. Moreover, the 
synchrotron emission is sufficiently weak at 150 GHz not to 
require any cleaning, which removes the need for low fre- 
quency channels to measure it. Experiments like EBEX and 
BICEP already match such conditions, not only because of 
the design choices, but also because their target areas inter- 
sects the PGMS strip (|Grainger et all 120081 ; IChiang et ail 
l2009h . 

A deeper detection limit, down to Sr = 2 x 10 , could 
be reached by a sub-orbital experiment observing the same 
area but with the CMB channel shifted to 70 GHz. The 
inclusion of a channel at a lower frequency would be required 
to measure the synchrotron component. Finally, detections 
limits down to Sr = 5 x 10~ 4 coiuld be achieved by observing 
at 70 GHz in a large area of 50° x 50° having the PGMS 
foreground levels. The location of the most suitable region 
must be determined, a task that can be accomplished by 
the forthcoming large foreground surveys like the S-band 
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Polarization All Sky Survey (S-PASS), or the C-band All 
Sky Survey (C-BASS). 

These limits are comparable to the goals of the space 
missions currently under study such as B-POL an d CMBPol 
|de Bernardis et al.1 I200II iBaumann et al.ll2008l ), but with 
the significant advantage that such an area is still sufficiently 
compact to be observable by a sub-orbital experiment. The 
limit r ~ f x 10 -3 is an important threshold for the in- 
flationary physics since it is about the lower limit of the 
important c lass of inflationar y models with low degree of 
fine tuning l|Bovle et alJliooil ). Our study shows that this 
threshold may be reached with an easier and cheaper sub- 
orbital experiment rather than a more complex space mis- 
sion, making this goal more realistically achievable with a 
smaller budget and in a shorter time than that required to 
develop space-borne equipment. 

The PGMS data will be made available at the site 
http://www. atnf. csiro. au/people/Ettore. Carretti/PGMS 
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